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THE  THEORY  OF  GAIN,  INTENSITY  AND  POWER 
OF  GAS  LASERS 

Kao  Cho#  and  Chao  Shu  Tao 

ABSTRACT 

The  particle  thermal  velocity  distribution  is  introduced 
into  the  speed  equation  to  obtain  the  relations  between  gain, 
intensity  and  power  of  gas  lasers.  This  paper  reached  the  same 
conclusions  as  reference  [1]  under  uniform  broadening  conditions. 
The  results  were  simplified  to  become  the  well  known  non-flowing 
gas  laser  equations  when  the  flow  speed  is  zero. 

In  the  analysis,  the  change  in  the  reflectivity  of  the 
mirror  and  the  degree  of  excitation  as  a  function  of  time  of 
gases  in  the  upstream  region  of  the  optical  chamber  were  consi¬ 
dered  with  respect  to  both  continuous  and  pulsed  gas-flowing  CO 2 
lasers . 

I.  INTRODUCTION 

In  the  analysis  of  a  gas-flowing  CO^  laser,  under  the  condi¬ 
tions  that  gain  is  equal  to  consumption,  the  solution  to  the  speed 
equations  put  on  line  [1,2]  applies  only  to  the  uniform  broaden¬ 
ing  situation.  Under  conditions  that  both  uniform  and  non-uniform 
broadenings  exist,  reference  [3]  made  an  analysis.  In  addition, 
the  semi-empirical  analyses  of  gas  flow  lasers  [4,5]  suggested 
that  It  was  possible  to  apply  some  theory  in  the  treatment  of  non¬ 
flow  gas  lasers  in  the  explanation  of  experimental  results  with 
gas  flow  lasers.  However,  theoretical  results  were  already  avail¬ 
able  [1-31  for  gas  flow  lasers.  But  they  cannot  be  reduced  to 
* 
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non-flow  gas  laser  equations  when  the  flow  speed  approaches  zero. 

The  Introduction  of  thermal  velocity  distribution  to  the 
speed  equation  will  enable  us  to  obtain  the  theoretical  relations 
between  the  gain,  intensity  and  power  of  gas  flow  lasers.  When 
the  flow  speed  equals  zero,  these  relations  reduce  to  the  familiar 
non-flow  gas  laser  equations.  Under  uniform  broadening,  the 
results  are  the  same  as  those  reported  in  references  [1,2].  It 
must  be  pointed  out  that,  vigorously  speaking,  particle  thermal 
velocity  distribution  cannot  be  brought  into  the  speed  equations. 
However,  this  technique  very  simply  provided  desirable  results 
and  will  serve  as  a  basis  for  further  analysis. 

II.  ASSUMPTIONS  AND  BASIC  EQUATIONS 

Let  us  assume  the  following;  The  optical  axis  is  perpendi¬ 
cular  to  the  direction  of  gas-flow,  gas  passage  is  rectangular 
in  cross-section  and  plane  parallel  mirrors  are  placed  on  either 
side  of  the  passage  as  shown  in  Figure  1.  The  variations  of  gas 
parameters  u,  p  and  T  Inside  the  optical  chamber  can  be  neglected. 
The  effect  of  the  boundary  layer  is  also  negligible  [7]  . 

The  pump  region  and  the  optical  chamber  are  separated  in 
order  to  analyze  the  "flow  broadening"  pulse  width  of  the  pulsed 
pump  upstream  of  the  optical  chamber. 

The  working  energy  levels  of  the  CC^-^  laser  system  can  be 
described  as  the  five-energy  level  model  shown  in  Figure  2.  Energy 
level  1  includes  the  symmetric  and  bending  levels  of  CO,,.  The 
number  of  particles  transported  from  energy  level  i  to  j ,  due  to 
inelastic  collisions,  in  a  unit  time  is  K. ,N, ;  N  is  the  density 
of  particles  at  energy  level  i.  Speed  (unit  sec  )  satisfies 

the  following; 


K»,  Kn,  K *  »  Kn  >  Kt;Q  -  1,  2,  3) 


(2.1) 
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With  constant  u,  p  and  T,  the  non-constant  speed  equations  and 
the  radiation  exchange  equation  are: 


Figure  1.  Schematic  diagram  Figure  2.  Energy  levels  of  the 
of  a  transverse  gas  flow  laser  CCU-^  laser  system  and  relaxa- 

and  coordinate  system.  tion  processes. _ 

(pump  zone'*. optical  chamber  -  o 

1 — gas  flow;  2 — pump  zone; 

3 — laser  direction 

where  t  is  time,  u  is  flow-speed,  w^(i=l,2,3)  is  pump  speed; 
c  is  speed  of  light,  h  is  the  Boltzmann  constant,  v  is  light 
frequency,  J  is  radiation  intensity;  f1  (1=1,2)  is  the  fraction 
of  particles  in  with  laser  action,  Bi1  is  the  speed  of  excited 
particles  with  Doppler  frequency  v*  under  the  influence  of  radia¬ 
tion  with  frequency  v. 
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(2.1*) 


g  (1=1^2)  is  the  statistical  weight  of  energy  level;  Av^ 
is  the  full  width  at  half  peak  of  the  uniform  broadening  curve. 

I  is  the  direction  of  light  propagation.  The  intensity  satisfies 
the  following  boundary  conditions 

y  —  0  j;  -  RJi  ,  R\  —  1  —  0\  —  *\  1  . 

y-L,  JC,-  R,Ji R,  -  1  ~  «.  -  *«  * 

Ri*  ai  and  ^i  2)  are  the  reflection,  absorption  and 

transmission  coefficients  of  the  mirrors  on  both  sides  respective 
ly.  J  and  J~  are  the  photon  energy  flows  in  the  positive  and 
negative  y  direction  and  we  have  J  —  J*  + 


III.  SOLUTION 

dj_ 

On  the  left  hand  of  (2.3)>  with  the  exception  of  ether 

terms  are  small  enough  to  be  negligible.  Therefore,  from  (2-3) 
and  (2.5),  we  get 

(L'  ±<.B»ftN,- B»UN{)<ly-=±\nRxR,  (3>1) 

L,  J*  e  2  Lt  v  J  ' 


Because  [1]  —  2>/R/(l  +  R~),  R  —  mm  (Rt,  Rt),  }  there 

fore,  the  variation  of  J  along  the  y  direction  can  be  neglected 
when  R^S0.6  .  From  (2.3)  and  (2.5),  we  obtain 


*/  -  f  r  -  (*'»/A.  -  B'dMJly  +  J;  -  K) 

L,i*  c  L , 


(3.2) 

(3.3) 


The  g  here  is  the  average  gain  coefficient  with  respect  to  y 
Equation  (3-1)  is  the  condition  required  for  gas  flow  laser  reson 
ance.  Generally,  and  g  vary  with  x. 
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The  solution  to  the  speed  equations  is  obtained  as  follows: 
First  take  the  integral  average  of  (2.1)  with  respect  to  y  and 
then  carry  out  these  mathematical  transformations: 


(2.1)  is  then  changed  into: 


X 

u 


(3.4) 


+  (*»  +  +  S,KW )  4-  S2K*Kunt 

Of1  of 


-  2  «/.  +  +  5c  [(K*  -  Ka>  If  4-  KUK  J 

"  Of  1  Of  J 

~i  +  (*„  +  5,**  4-  S,KW)  4-  - 

Of1  Of  ... 

-  5,  K*  S  ^  +  5.*»  ( +  fiutg  ) 

,  i  - 1  Of  'Of  / 

where 

»*  ”  "l  +  »j,  »j  “  5,n#  4-  S„g ,  «i  —  S2n»  —  S«f 


(3.5) 


5. 


-,5> 


fln/i 


(3.6) 

(3.7) 


Bi i/a  4-  B{j/i  ’  flji/j  +  Sn/i 

..  -  f  (*•  -v, 

L,  J» 

Let  e  -  *(?)*'*»  and  then  obtain  the  solution  to  (3.5)  as: 


»*  -  «»,  +  5.*»»«  4-  2  rZT  ((«'»  +  w>  ~  + 

/  -  1 

4*  (jl,  —  SiK*  —  5aKip)nS(0  +  So*[(Kio  —  K«)  —  ih  +  $)»*#]*  #i) 
».  —  nJr  +  S,  gnit  4-  2  ~.  f  ]  K«'>  “  **»>»)  +  “  Ksi)*»(0 

f«u) 

4-  S,M(C)  +  Sog[K>'  -  (*-  + 


(3.8) 


The  superscript  0  on  the  right  represents  the  distribution 
of  the  corresponding  quantity  at  g  =  0  (i.e.,  x=0).  6  is  a  con¬ 

stant  and  w^  is  another  constant. 

A,  ,  -  1  {(K.J  4-  5,K„  +  S>K,a)±V (/'Ji  4-  S^u  4-  Sj^ic)1  —  «S,KUA:U  ] 
2 


5 


»»« 


nu 


_ (^m  Ki\)8  +  K»K\t 

S'  +  (Ku  +  SxKa  +  S1Klt')8  +  StKuK*  * 

_ *»(*».  +  S') _ 

S !  +  (K«  +  S\Ki*  +  SjK^S  "+”  SiKnKii 


n>. 


Ly„ 

<iS,Ku£i  * 


o<£<£> 
£  >£> 


i>0 


Sjkiau,i  "+■  ^  \  wt 

_ <  - 1 

S,KaKv 

\  0 


0<£<£, 

£>£> 


(3-8)  is  the  solution  when  A^  A Similarly,  a  solution  when 
A  =  A2  can  be  obtained.  It  will  not  be  discussed  here.  The 
initial  distribution- at  x  =  0  is  explained  as  follows:  (1)  ini¬ 
tial  distribution  can  be  obtained  from  the  solution  to  (2.2) 
under  radiationless  condition  Ji  0.  The  radiationless  solution 
is  not  going  to  be  written  here.  Under  continuous  pumping  at  the 
upstream  of  the  optical  chamber,  i»?(C)|t-o—»?(<)  |»^—  constant.  When 
pulsed  pumping  is  used,  »!(*)!«-•  varies  with  t.  (2)  the  radia¬ 
tionless  solution  usually  generally  does  not  satisfy  (3-1).  How¬ 
ever,  since  — ; — 2L L  ■-  «  Kil,  -^V  ,  the  inelastic  exchange 

r(l  —  RiRi)  Wi  ^ 

of  collision  within  the  time  period  f(i— R  ^  cannot  be  completed. 
In  addition,  since  — ; — — -  «  L„  at  £  *  0  we  get: 


n\  -  «KO  +  »?( O,  n\  - 

where  ”?(f)  “  "?(r)  "  l'*  2»  3) 

radiationless  solution.  (3)  g 


;herefore 


Ku*5(f)  ™  ^d°jCO. 


•  "SCO 

which  can  be 
enerally,  since 


(3.9) 

obtained  from  the 

,  K\t  Kn,  — , 


From  (2.2)  and  (3-6),  gl  can  be  obtained 


Si  . 

hv 


SjWi  *4*  j  //2,  —  StfCd  “t  &21  t  ^2^23  t 


_2 _KtL  ^  f  r2(v  —  f0)v  In 2 1  a  1 
itAvh  hv  exP|  j  Avg  J  J 


A  L  Avs  J  J  hv 


(3.10) 
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where 


«Av. 


/  -  4-  S2Ka  -  S}Ku  (3.10), 

— —  ="(«'»  +  wv  —  ~l~2)  2  “'<»)  2  - T  — 

hv  '  i-1  '  *, 


>-IJ) 


—  5|/(«'io  4*  «/»  —  1/fl*,)  4"  l$l^.'l(l*  —  K„)  —  Sl/Kjj]n5o(f) 

4-  iS,S,KaKu  -  SJCX,  -  S,KU  -  SjK.JKXO) 

L  „  -£ll  {(Kj,  4-  S2K„  +  S\Km  4 -  S  -  £,£„*„  4-  S]fnt,') 

hv  2  Bufj  l 

-  2  (S,K.2[Kv  -  (X,  +  S)n)t]  -  SjUKu,  ~  Ka) 

I  •  1  (;^i .  **f  *•* 

#•14) 


-  (A,  4-  *)»»,])} 


In  the  derivation  of  (3-10),  it  has  been  assumed  that  a 
quasi-equilibrium  was  reached  and  partial  Maxwell  velocity  dis¬ 
tribution  was  established,  i.e.. 


u>t  —  u/ exp 


n,  —  n  it  exp 


where  vq  is  the  center  frequency  of  the  Doppler  line  and  Av^  is 
the  full  width  at  half  peak  in  the  Doppler  line. 


IV.  GAIN,  INTENSITY  AND  POWER 


Using  (3*10)  under  conditions  that  the  light  frequency  equals 
the  Doppler  center  frequency  (i.e.,  v  =  vq ) ,  we  get 


e 


*.('  exp(-=V)  , 
*  J— 1  4-  r2  4-  Hi,  Z 


(4.1) 


_  Avy  / 1  1  2(v‘  —  ve)  , 

7  ~r~a~  V  In  2 ,  *  —  - *S  ext 

where  Av0  Avw  is 

When  non-uniform  broadening  dominates. 


the  probability  integral, 
i.e.,  n  -*■  0,  (4.1)  becomes 


r  -  *.(?>  o 
Vi  +  7777(f) 


(4.2) 
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When  uniform  broadening  dominates,  i.e.,  n-*-00  (4.1)  becomes 


t 


slLSL- 
i  +  ///,(*)  * 


(4.3) 


It  is  defined  as  the  partial  saturation  intensity  Kq  and  gQ 
are  non-uniform  and  uniform  broadening  saturation  gain  coeffi¬ 
cients,  respectively.  When  the  gain  is  equal  to  the  loss,  the 
above  equation  also  applies.  It  becomes  the  familiar  theory  [6] 
on  non-flow  gas  lasers.  When  l,S  1 

1  _3  ttAvn  chvS,  KiiKuKio  +  +  (AT,,  +  Kj\  +  Kii)^nl 

2  B1Jl  SiKuKu  +  ( Ka  +  SxKa  +  SaK»)* 


Equation  (4.4)  omits  6  as  a  higher  order  term.  From  (4,1)- 
(4.3),  it  is  easy  to  obtain  the  radiation  intensity  in  the  opti¬ 
cal  chamber.  The  transmitted  intensity  Jt  is 


]. 


Jiti  + 


glL}(ttVRt  +  tx*/ R  i) 

cVrx  +  7^)0  -  Vr#,) 


(4.5) 


Integrating  J  with  reppect  to  x  and  t,  we  get  the  power. 

For  single  end  output  with  the  other  end  as  a  no-loss  total  reflec¬ 
tion  mirror  (i.e.,  R2  =  1)>  we  get 

v_ . .  + %_.,(■  r^)W|  c«) 

where  V  ■  is  the  volume  of  the  optical  chamber.  For 

ij  »  1  ,  (4.5)  becomes 


1  77F  ^r"V  T  ft  I  +  \"  I  1 

L,(C  —  C«)(«i  +  *»)  J*  I  L j  Jc.  J*  I  h) 


where 


-rr~  \  +  laRv 


(47) 


IT - - - l,'*td'd c| 

Ll  L  I'  I «_ 

l-iywtL 

.  JC,  j* 


(4.8) 
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I*  and  g*  can  be  considered  as  the  gain  coefficients  of  uniform 
so 

saturation  intensity  and  non-uniform  saturation  intensity,  res¬ 
pectively.  They  are  average  values  with  respect  to  time  and 
space.  (4.7)  is  in  the  same  form  as  the  Rigrod  equation  for  non¬ 
flow  gas  lasers.  Substituting  (3. 10)  into  (4.8),  I*  and  g*  can  be 

s  o 

obtained : 

IT  -  -h-Ufl  ( (K*  +  SjKa  +  S,KU  +  8  -  S}Kunu  4- 

2  BJjL ,  l  .  8 

^  ^  ~~p~ - r Cl*  +  <*)">,] 

i  -  s(M 

-  s,/[(k„  -  -  a,  +  o-9> 

The  expression  for  g*  is  going  to  be  omitted.  The  t*  corres¬ 
ponding  to  the  maximum  power  out  is  obtained  based  on  dP  _Q 
For  6  =  0,  we  get  9/* 

~  -  ~  '*  IgTL,  +  lo(l  >-  -  /,»)]  (4  10) 

“ 1  a  1  +  *7 


From  (4.7)  and  (4.10),  the  maximum  power  output  P*  is 


P * 


virory 


on1 


Ljfli(l  a,  «i(l  at  —  1 1)  L ilj(i 


"r-^TT  f  ^  I  010 

^  —  Cc)  Jc, 


The  above  equation  has  the  same  form  as  the  Rigrod  equation 
[6].  For  7  "K  1  (non-uniform  broadening),  the  same  discussion 
can  be  made. 


V.  ANALYSIS  OF  PARAMETERS 

2 

I  is  proportional  to  the  square  of  the  pressure  P  .  It  is 

not  related  to  the  pumping  rate  and  the  degree  of  excitation  of 

the  gas  at  the  inlet  of  the  optical  chamber.  From  (3-10)  we  know 

that  I  decreases  monotonically  with  —  from  its  maximum  I 
s  us 
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max  at  x  =  0  to  I  (see  equation  (4.4)). 


—  /,1,-B  =  (*»  +  +  S.K-  +  *)  5r'^rL 


Oi) 


Increasing  the  flow  speed  (leaving  other  parameters  un¬ 
changed),  we  found  that  I  remained  the  same  but  I  and  I 

s  y max  s  s 

increased  significantly  (see  example).  The  variation  in  reflect¬ 
ivity  R1  =  (R-^° )  ^  significantly  affects  Ig.  The  reasonable 
limits  for  6  was  derived  to  be: 

\S\  <mina„  X,)for|»|  -  0(K*)  (5.2) 

K  K 

Prom  (4.4)  we  know  that  e  **  0  •  This  is 

because  when  $<o  R^  increases  monotonically  with  —  and  t^ 
monotonically  decreases.  When  t^  is  reduced  to  zero,  light  radia¬ 
tion  is  forced  to  stop;  when  5>0  >  or»  )//!*-*. 

In  gas-flow  lasers,  K  and  g  are  parameters  of  I  which  are  re- 

o  o  s 

latea  to  the  mirror  surface  conditions.  Only  when  u  =  0  and 
6=0,  these  parameters  will  have  the  physical  meaning  of  corres¬ 
ponding  parameters  for  non-flow  steady-state  gas  lasers.  It 
demonstrates  the  characteristics  of  the  non-steady  state  nature 

of  gas-flow  lasers.  The  parametric  dependence  of  I*  is  identical 

s 

to  that  of  I  . 

s 


If  the  pumping  rate  w^  is  proportional  to  p,  then  the  non¬ 
saturation  gain  coefficient  K  is  not  related  to  p  and  g  is 

o  o 

inversely  proportional  to  p.  Kq ,  gQ  and  w^  are  proportional  to 
the  degree  of  excitation  of  the  gas  at  the  inlet  of  the  optical 
chamber.  From  (3-10)  and  (4.4),  we  know  that  for  pumping  upstream 
of  the  optical  chamber,  for  6  >  0,  gQ  decreases  monotonically  and 
for  6  <  0,  gQ  increases  monotonically.  For  pumping  inside  the 
optical  chamber,  we  get  from  (4.4) 

K#  ”  («'»  +  w*—  2  **'«) 

.  By fdSiKuK,'  -4-  (K»  +  S,yp  +  _  (5  3) 

cSAK„KuKlt  4-  l K„KU  +  (Xu  +  Ku  4-  K^K^S) 
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Therefore,  for  the  same  pumping  conditions  ,  £»!«><> <(pr  « 

Power  P  is  proportional  to  p,  u,  pumping  speed  and  the 

degree  of  excitation  at  the  inlet  of  the  optical  chamber.  When 

reflectivity  is  constant  and  6=0,  the  maximum  power  output  P* 

2 

is  proportional  to  p  and  mirror  area. 

VI.  EXAMPLE 

For  a  C02/N2/He  gas  mixture,  the  speed  constants  are  tabu¬ 
lated  as  follows  [8] : 


speed  constants  at  T  *  300°K 


Uu  /N 

He 

k23/p 

K32/K23 

(T0RR-1sec~1 ) 

0 

8.8  x  102 

1.23  X  102 

1.67  x  104 

j 

0.3 

0.5 

1.2  x  103 

1.9  x  103 

1.03  X  102 

9.6  x  101 

1.16  x  104 

8.3  x  104 

;nco2/ns2 

From  information  in  reference  [8],  we  get 

Si  Afo  2(2/  +  Q 
S<,  T 

S^/S2  is  also  a  function  of  temperature  and  degree  of  excitation. 
When  T  =  300°K  and  degree  of  excitation  is  not  too  high,  ^  0.(M 

The  variation  of  I  with  ^  <  ^  is  shown  in  Figure  3. 

s  3  1  2 

I  approximately  is  in  the  10J  -  10  watt/cm  region  which  is  con- 

*  al 

slstent  with  experimental  results^.  The  calculated  curve  of  Ic 
proved  the  analytical  conclusion  given  in  the  previous  section. 
When  8  —  lj  ,  I  is  10  times  larger  than  I,!*-* 
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Figure  3-  Variation  of  I 
with  ii  s 

M 

r-  30UK,  p  -  30  tnjp,/s./H< 


1/VS 


Figures  4  and  5  show  the  change  of  g  and  I  with  and 

ut  O  u 

y  ,  respectively.  The  pulsing  time  of  the  pump  vs.  time 

*->  o 

can  be  expressed  as  the  variation  of  with  —  .  gQ  and  I 

are  like  "waves"  propagating  in  the  right  direction.  Straight 
lines  —  and  =  1  are  the  front  and  back  surfaces  of 

the  wave  where  b  —  — .  .  From  Figures  4  and  5,  we  know  that 
along  a  straight  line'  y  ~  b^-  ~  H  —const  o<  h<\,  the  varia¬ 
tion  of  gQ  is  slower  than  that  of  I. 

For  continuous  pumping,  when  R^  =  const,  go*  monotonically 
gradually  decreases  with  ih*.  .  When-^  —  10>  gf  is  approximate¬ 
ly  0.8.  When  $  =» -i*  g  *  more  significantly  monotonically  declines. 

.  5  O  —X  _ 

When  —  —  10  «*  is  about  0.1.  When  8  *■  — y L 1  ,  g  *  increases  with 

J  „  M  *  5  Q 

— .  and  I  *  decreases  with  *»£-  .  Note  that  g  *  and  I  *  do 

u  S  u  dO  S 

not  have  any  obvious  physical  meanings.  From  Figures  4  and  6,  it 
can  be  seen  that  when  R  is  constant,  the  variations  of  gQ  and  gQ* 
with  ^  are  slow  regardless  of  whether  continuous  or  pulsed  pumping 
is  used.  This  serves  as  a  theoretical  basis  for  the  calculation 
of  gas  flow  laser  power  using  an  average  value  of  the  gain  coeffi¬ 
cient  and  the  non-flow  gas  laser  equation  as  suggested  in  refer¬ 
ences  [5,73- 

Gas  laser  power  calculated  here  agrees  with  that  reported 
in  reference  [1]  (see  Figure  7).  represents  the  fraction  of 

effective  vibrational  energy  taken.  For  the  calculation  para¬ 
meters  in  Figure  7,  when  «,  5  ,  the  effective  energy  is  com- 

U 

pletely  used. 
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m 


Figure  6.  Unified  average 
non-staturation  gain  g  *  vs. 
A2x  ° 

u 

T=300°K  p=30  torr 

C02/N2/He-i/4/5 


Figure  7. 
power  g_  vs 
uA 


Gas  laser 
.  A  zx 
u 


2.  6=0,  3 


output 


T=  3 0 0 ° K ,  p=30  torr  . 

CO,/N-/He-l/4/5  7^jr, 

SX=0  *  2x1°:3=">'1  "2  ■ 

particle/cnr 
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CONTACT  PROBLEMS  OF  LONG  RIGID 
FRAME  FOOTING  ON  ELASTIC  FOUNDATION 

Zeng  Xin-Chuan* 

(The  Seismological  Brigade  of  Wuhan,  National 

Seismological  Bureau) 

In  reference  [1],  Galin  solved  the  contact  problems  between 
two  flat  bottom  pressure  heads  under  center  load  and  elastic 
semi-flat  foundation  with  friction  present.  The  distribution  of 
stress  along  the  contact  surface  was  obtained.  In  this  paper, 
the  author  used  the  Muskhelisbvili  method  [2]  to  attack  the  con¬ 
tact  problems  of  long  rigid  frame  footing  under  off-center  load 
on  an  elastic  foundation.  For  simplicity  in  the  calculation,  it 
was  assumed  that  the  contact  between  the  footing  and  the  founda¬ 
tion  is  frictionless. 

The  action  of  an  off-center  load  (Figure  1)  is  equivalent 
to  the  combination  of  the  actions  of  a  center  load  and  a  force 
pair  (Figure  2).  The  boundary  conditions  of  center  loading  are: 

f —  0  y  —  0  lx)  <  00 

<r>-0  y-0  Ixj  <  a  \x\  >*.  (D 

'‘'■'f  (IS)  y  — 0  a  <  jxj  <  b 


Figure  1  Figure  2 

If  ~ 

This  paper  was  received  on  November  7>  1978. 


15 


For  the  force  pair,  the  boundary  conditions  are: 


x*y  “  0  y  —  o  |*J  <  oo 
<r,-0  y-0  |*|  <*  |*j  >b 

»  —  jx—dx  y  —  0  «<|*|<A 


(2) 


where  0  =  d/b  is  the  rotation  angle  of  the  footing. 


The  displacement  and  stress  of  the  foundation  can  be  expressed 
by  the  complex  function  <j>( z): 


2/*(«  +  «V)  —  *?(*)  +  <p(»)  +  ( =  —  *)<p'6f) 

—  <r,  +  2 it,,  —  —  2{(*  —  i)<p''(*3  —  <?'(*)  ~  pT*!) 

<r,  —  it,,  —  <p'G?)  —  <?'(*)  +  C*  —  *)<?>"(.*) 


(3) 


The  solution  of  <j)(z)  which  satisfies  the  boundary  condition 
(1)  is: 


<Jj(*)  —  ‘P~  — 

v  2**00 


o> 


where  *  “  »  X(.*)  **•  C*1  —  *1')l/,7 


.  The  solution  of  <f>(z) 
which  satisfies  the  boundary  conditions  (2)  is 

i  a:  (  ** - T  +  ^*)  1 

- kr— }  C5) 


Under  center  loading  conditions,  the  displacement  of  the 
footing  can  be  obtained  from  (4)  and  (3): 


terVti,  {ogA  a<  1*1  <*  (6) 


where 


A  -  (b*  -  f-arctg 


b  si 


sin  (0 


( \b *  co*1  to  —  aiy/l  ’ 


u  ™  arc  aa 


(*’  -  «*) 


l/» 


For  footing  under  the  action  of  a  force  pair  M,  the  angle  of 
rotation  is: 


-  (1  +  QA/ 

*2*m<*j  “ 


C7) 
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The  stress  distribution  on  the  bottom  of  the  footing  under 
center  load  conditions  can  be  obtained  from  (4)  and  (3)  as: 


r(x2  -  s’yxv  -  x1)* 


I  r**  “  0 


a  <  |*|  <  b  | 


Under  the  action  of  the  force  pair,  from  (5)  and  (3)>  we  get: 


<ryU“B  -  *2)  (xl  -  -x'Y*  '  "  "  *  "  '  f  W 

As  for  the  stress  distribution  of  the  foundation  under  center 
loading,  we  obtain  the  following  by  substituting  (4)  into  (3): 
W*W)-*  [c-  A(x,  y)  +  l)y  CO.  -?i— .?!  ±J»L  | 


X2  -  y  („2  +  A2) 


«  <  |x|  <  b 


+  (y£(x,  y-)  ~  x)sin 


Po  +  <jpq  +  <p!  +  go? 


y  (p^p2p•lpD_,',  [(^(x,  y)  +  l)y 


P»  +  <Po  +  <Pl  +  <Pl 


—  (y5(x,  y)  +  x)  sin 


Po  <Po  <pj 


r,,  —  —  ^  CpJp6p!pD-1/j  [^(x,  y)  sin  -g--— 

+  B Or,  y)  cos  ?’■] 


n*  4>  4-  m! 


where 


[(x  -  <»)J  +  y2]1'2  PC  -  [(x  +  a)*  +  y’l1'2 


<pi  —  arc  tg 


qPo  —  arc  tg 


aJ  t-  y  i  \ 

x  +  a  J 


Pi  -  ICx  -  by  +  y2]1*  p?  -  [(x  +  +  y2]1 


<pl  »•  arc  tg 


*  —  b 


tp]  —  arc  tg 


x  +  * 


/<(x,  y)  -  [(x2  -  y2  -  *2)(x2  -  y2  -  „2)  -  4x2y2]((x2  -  y2)2  -  4x2y2  -  *2A2] 
+  8x2y2(x2  -  y2)[2(x2  -  y2)  -  (o2  +  A2))/[(x2  +  y2)2 
-  2(x2  -  y2)*2  H-  A‘][(x2  +  y2)2  -  2(x2  -  y2)«2  +  o'] 

0(x,y)  —  2xy[(«2  +  £2)(x2  +  y2)2  —  4«2£2(x2  —  y2)  +  «V(«2  ■+•  £2)]/ 

C(x2  +  y2)2  -  2(x2  -  y’W  +  *«j[(x2  +  y2)2  -  2(x2  -  y2)*2  +  o'] 
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Under  the  action  of  the  force  pair,  we  get  the  following  by  sub 
stituting  (5)  into  (3): 


p*  +  vi  +  <p! 


CpWpW)-"  {[(**  -  -W*.,)  -  4*]y  cos  ?9  ±  S&±  SL  ± 
+  [2(x1-  /-!(«■  +  i1))  -  (i'  -  W(i.y)  |  aa2l±2!±£!±S 

-  -  T7r^C^W)^{[C^--0,c(x,y)+4,]yco«  ^ ±jeL±jeL± 

—  a  )  1  2 

—  — /  —  y(.’  +  i')j  +  (O'  —  jOVDCt,  >)  j  da  Tl  +  2l± Si +  '• 

r„  _  [c<,(  ,)„i.  Ti-1-.tj+f.'f  Ti. 


i  +  +  j>!  +  <pi 


+  D(.x,  y )  cos 


Po  +  <p{  4-  qpf  | 


where 


p1±jpL 


-  *K*a  -  yQ1  -  4/  -  W  +  /Xo*  +  **)  + 

[(*•  -y  -  ay  +  4x,y‘]K**  -  y*  -  #*)*  +■  4xyj 

y[4x«  -  (xJ  -  y1)1  -  (*>  +  y'X*1  +  **)  -  -V] 

y)  -  t(xi  _/  -  b*y  +  4*yuc»*' -  7  -  77  + 


_n  (14) 
JEll 


Figures  3  and  4  show  the  stress  distribution  a ^  on  a  few 
horizontal  surfaces  under  center  loading  and  force  pair  condi¬ 
tions  respectively.  Each  figure  only  shows  half  of  the  entire 
picture.  The  other  half  of  Figure  3  is  symmetric  to  what  is 
shown.  For  Figure  4,  the  half  nut  shown  is  asymmetric  to  what 
is  presented. 


Figure  3 


Figure  4 


From  (6),  (7)  and  (8),  we  can  see  that  when  a  =  0,  the 
following  holds  for  long  solid  rigid  frame  footings: 


e 


(*  —  Opto 


» 


Cl  ±  OAf 
Ixfib1  ' 


v 


y^ 


4xfi 

_ P 

*(b>  -  t2)1'*’ 


1*1  <  b, 

|*|  <b. 


Our  results  are  consistent  with  those  found  by  Galin  [1],  Musk- 
helishvili  [2]  and  Sneddon  [ 3 ] - 


From  (9)  and  Figure  4,  we  see  that  changes  signs 

between  -b  and  -a  and  a  and  b.  This  implies  that  tensile  stress 
will  appear  at  some  location  on  the  bottom  of  the  footing.  When 
the  footing  does  not  adhere  to  the  foundation  or  when  the  adhe¬ 
sion  is  less  than  the  tensile  force,  the  footing  will  be  detached 
from  the  foundation.  Therefore,  in  order  to  avoid  such  a  detach¬ 
ment  effect,  the  load  and  its  corresponding  force  pair  must  satisfy 
certain  relations. 


From  (8)  and  (9),  we  get  the  stress  distribution  on  the  bot 
tom  of  the  footing  under  off-center  load: 


_ 2£ _ 

,<y  -  aO^O2  -  *2)’/> 


2M 

*(A2  -  ,2) 


**  -  y  (V  +  *2) 

C*2  -  *2)‘/2C  b1  -  x1)1'1  * 


a  <  1  *  |  <b 


It  is  well  known  that  in  order  not  to  allow  detachment  of  the 
footing  from  the  foundation,  it  is  necessary  to  make  the  normal 
stress  on  the  bottom  of  the  footing  be  pressure  stress.  Therefore, 

2M  *a-y<y  +  62)  _ 21 - ,  a  <  1*1  <  * 

(A2  -  a2)  (*2  -  -WA2  -  *2)‘/>  (*2  -  «2)1/1C^2  -  *2)‘/2 

and  M  <?(*  +  *)/  2 
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If  the  action  point  of  load  is  a  distance  r  off-center  (see 
Figure  1),  then  M  =  rp.  From  this  the  equation  that  r  must  satis¬ 
fy  so  that  the  footing  and  the  foundation  will  not  separate  is: 

r  +  £)/2  (17) 

For  long  solid  frame  footing,  a  =  0,  and  the  above  becomes: 

M<pb/ 2,  r<b/2 

This  is  consistent  with  result  obtained  by  Muskhelishvili 
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